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Abstract 

A detailed comparison of the expressions for the decay widths obtained within the semiclassical 
WKB approximation using different approaches to the tunnehng problem is performed. The dif- 
ferences between the available improved formulae for tunneling near the top and the bottom of the 
barrier are investigated. Though the simple WKB method gives the right order of magnitude of 
the decay widths, a small number of parameters are often fitted. The need to perform the fitting 
procedure remaining consistently within the WKB framework is emphasized in the context of the 
fission model based calculations. Calculations for the decay widths of some recently found super 
heavy nuclei using microscopic alpha-nucleus potentials are presented to demonstrate the impor- 
tance of a consistent WKB calculation. The half-lives are found to be sensitive to the density 
dependence of the nucleon-nucleon interaction and the implementation of the Bohr-Sommerfeld 
quantization condition inherent in the WKB approach. 

PACS numbers: 03.65.Sq, 21.10.Tg, 23.60.+e 
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I. INTRODUCTION 

The Wentzel-Kramers-Brillouin or the WKB approximation , sometimes also known 
as the BWK j3| , the semiclassical approximation or the phase integral method has been 
widely used in the evaluation of the half-lives of radioactive nuclei. It was evident from the 
historical papers of Gamow [6] and Condon and Gurney 7(] that one could treat the alpha 
decay of nuclei in terms of the tunneling of a preformed a-particle confined to the interior 
of the nucleus, through the Coulomb potential barrier of the alpha-nucleus system. The 
WKB approximation which is really applicable when a problem can be reduced to a one- 
dimensional one was found suitable to evaluate the barrier penetration probabilities and the 
decay width in general was defined as a product of the frequency of collisions of the a with the 
barrier (the so-called assault frequency) and the penetration probability. The objective of the 
present work is to critically examine the decay widths obtained within some entirely different 
approaches to the tunneling problem, however, all working within a WKB framework. To be 

n 

specific, we examine four approaches; the first a two potential approach (TPA) [8], another 
a path integral method with Jost functions j^j, a third one using comparison equations to 



obtain improvec 
model (SAFM) 



WKB formulae and a fourth one involving a super asymmetric fission 



10[ |. It is gratifying to know that all the four approaches indeed lead to the 



same formulae for the WKB decay widths. A WKB expression for the vibrational energy, 
E^, emerges from the above comparisons. For the first three approaches, the vibrational 
energy (and hence the assault frequency) can be simply evaluated from the potential and 
the tunneling particle energy. In the SAFM however, the equation for E^, turns out to 
be a transcendental equation. However, this equation is often neglected by the SAFM 
calculations in literature and a fit to the vibrational energies is performed, leaving the 
calculation somewhat incompatible with the WKB framework. 

Besides this, the widths calculated within the SAFM often neglect the Bohr-Sommerfeld 
condition and the Langer modification which are essential ingredients of the WKB frame- 
work. The Langer modification is a necessary transformation while going from the one 
dimensional problem with x ranging from — oo oo to the radial one-dimensional tunnel- 
ing problem with r ranging from ^ oo. To demonstrate the importance of performing a 
completely consistent WKB calculation of widths, we perform a realistic calculation of the 
alpha decay widths of super heavy nuclei which are a topic of current interest. 



2 



With most radioactive decays occurring away from the extremes of the potential barrier, 
the standard WKB is found to be a reasonable approximation for such calculations. However, 
for specific cases, where the tunneling can occur near the top or the bottom of the barrier, 
the validity of this approach becomes questionable. Some attempts at obtaining improved 
formulae which can be used near the extremes of the barrier do exist {4, 9] and will be 
investigated by applying to realistic examples in the present work. 

The alpha - daughter nucleus interaction in the present work is described using folding 
model potentials with a realistic nucleon-nucleon (NN) interaction Such a model has 
been shown to be quite successful in predicting half-lives of unstable nuclei 12| . We perform 
calculations with the added ingredient of a density dependent NN interaction and find the 
results and fitted parameters sensitive to this input. In the next section, we perform a 
comparison of the different WKB approaches as mentioned above. In Section III, we briefly 
present the relevant formulae of the potentials used. The results are discussed in Section 
IV. 



II. THE QUASICLASSICAL ALPHA TUNNELING PROBLEM 

Starting with the description of a radioactive nucleus as a cluster of its daughter nucleus 
and an alpha particle, it is by now standard practice to study the alpha decay of nuclei as 
a tunneling of the a through the potential barrier of the alpha-daughter nucleus system. 
Though most quasiclassical approaches agree on the proportionality of the decay width to 
an exponential factor, namely, F oc e^^*^, where G is the famous Gamow factor, the details of 
the calculations vary depending on the approach used. Typically, one considers the tunneling 

V(r) 




FIG. 1: Typical potential for an alpha- nucleus tunneling problem. 
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of the a through an r-space potential of the form, 

V{r) = + V;(r) + + y^^' , (1) 

where Vn{r) and Vc{r) are the nuclear and Coulomb parts of the a-nucleus (daughter) po- 
tential, r the distance between the centres of mass of the daughter nucleus and alpha and 
u their reduced mass. The last term represents the Langer modified centrifugal barrier 

n 

PJ]. With the WKB being valid for one-dimensional problems, the above modification from 
Z(Z + 1) (/ + 1/2)^ is essential to ensure the correct behaviour of the WKB scattered 
radial wave function near the origin as well as the validity of the connection formulas used 



14| . Another requisite for the correct use of the WKB method is the Bohr-Sommerfeld 



quantization condition, which for an alpha with energy E is given as. 



r-2 



K{r)dr = {n + 1/2)tt (2) 

Jri 

where K{r) = ^^{\V{r) - E)\, n is the number of nodes of the quasibound wave function 
of a-nucleus relative motion and ri and r2 which are solutions of V{r) = E, are the classical 
turning points (as shown in Fig. 1). We shall now examine the decay widths obtained within 
different WKB based approaches to the above problem. 

A. The two potential approach 

n 

Starting with a typical potential as in Fig. 1, the authors in [8] consider the tunneling 
problem of a metastable (quasistationary) state and obtain a perturbative expansion for the 
decay width of the metastable state. The potential is split into two parts, V{r) = U{r) + 
W{r), where the authors first consider an unperturbed bound wave function ^o{r), which is 
an eigenstate of the Hamiltonian, Hq = — (/i^/2yu)V^ + U{r). W{r) is a perturbation and 
when it is switched on, $o('") is not an eigenfunction of H = Hq + W{r), but one rather 
has a wave packet which is expressed as an expansion in terms of $o(^) and the continuum 
wave functions, $fc('")- Once the perturbative expansion for the width is obtained in terms 
of the wave functions, retaining the first term and expressing the wave function using the 
semiclassical WKB approximation, the width is given as. 
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Ttpa{E) = — 



'"2 dr 



2/i [Jri k{r) 



-2 p k{r)dr 

e •^'■2 (3) 



where, k{r) = d^l'") — E\). The factor in front of the exponential arises from the 
normahzation of the bound state wave function in the region between ri and r2. Indeed, 
this factor is related to the so-called 'assault frequency' of the tunneling particle at the 
barrier. Expressing the time interval. At, for the particle traversing a distance, Ar as, 

v[r) nk[r) 

the assault frequency v can be written as the inverse of the time required to traverse the 



distance back and forth between the turning points ri and r2 as 



V 



h 

2/i 



'1-2 



dr 



'^{\V{r)-E\) 



-1 



(5) 



Thus, TxpAiE) = hue'"^^, where, W = J^^ k{r)dr. It is interesting to note that the 
above definition of the period T (which is twice the time spent between the turning points 
ri and r2), is formally similar to that of "traversal time in tunneling" as defined by Biittiker 
and Landauer [l^. They defined the time for a particle traversing the barrier and hence 
obtained a similar definition as that for T/2 above, but within the limits r2 to r^. 



B. Improved WKB widths 

The width obtained in the two potential approach is in accord with the widths obtained in 
4] and Qj for energies well away from the top or the bottom of the barrier. For example, in 
9j, considering a double humped (DH) barrier in one dimension (with cartesian coordinates) 

and using the method of comparison equations, the decay width for all energies except near 

the very top or the very bottom of the barrier was found to be, 

Tdh{E) = 2huDHH^ + e-2W] ^ 2hi/DHe-2W (6) 

where, udh is the assault frequency in the case of a double humped barrier. Replacing 
J^DH = for small values of the exponent, Tdh{E) indeed reduces exactly to Ttpa{E). 

1. At the base of the harrier 

The width for energies of the tunneling particle near the bottom of the well was found in 
9] to be, 

T%{Er,) = 2huDHH^ + a-'e-2W] ~ a-^T^n{^^) , (7) 



where, a ^ = (1/"^!) -\/27r [(n + 1/2) /e]""*"-^/^. The relation obtained in 8| using the TPA is 
somewhat different, namely, 

V2 



-plow 



TPAi^n 



TT 



(8) 



It is mentioned in 
time path integral methods [1 
expression for the widths in [l6 



16| obtained using complex- 



Bfl that the above result coincides with that in 

The prefactors in front of the exponential factors in the 
17| depend on characteristic parameters of instanton tra- 
lectories, however, the above modification in ([8]) does agree in form with that obtained in 
16 1 . Though the TPA approach and that of comparison equations in jo] do give the same 
expressions for the widths away from the extremes of the barrier (TTPAiE) = ToniE) was 
seen above), they do not seem to agree on the situation at low energies. Eqs ((Tj) and ([8]) 
seem to give very different results. For example, for n = 0, 1 respectively, one finds that 
r^^(E) = 1.075TDHiE) and T^EmiE) = 1.027TnH{E), whereas, P^PAiE) = AASATrpAiE) 
and Tiji^AiE) = 7.2A7Ttpa{E) respectively. In [18|, one can find yet another expression for 
tunneling near the bottom of the barrier such that the tunneling probability at the base of 
the barrier vanishes exactly. 



2. At the top of the barrier 

In js], the decay width within the WKB approximation at the top of the barrier was also 
evaluated and found to be, 



Y^top 
DH 



(E) 



dE 



(9) 



where, Tdh is the period for moving back and forth between the two humps and = 
argT{^ - iW/tt) + f [ln(W/7r) - 1]. A very similar formula was also found in , however, 
with a difference of a sign in the denominator. Having derived the expressions for the Jost 
function of a radial barrier transmission problem (with one hump) by the path integral 
method, the authors obtain the following expression for evaluating the width at the top of 
the barrier: 

2h[{l + e-2W^)V4 _ (1 + e-2^)-V4] ^^^^ 



Froman \ ) 



T 
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The numerators in (fTOj) and (Q for energies near the top of the barrier are almost equal, 
i.e., 2n[{l + e-^^y/^ - (1 + e-^^)"!/^] ^ mn[l + e'^W]. With the period, Tdh being 
twice as much as T and the factor = —2 a, Eqs and (ITU]) indeed agree but up to a 
sign in the denominator. In Eq. (5.7) of ^|, however, one can see that there exists a choice 
for the sign appearing in front of da/dE and the authors choose the negative sign without 
any particular justification. We shall later notice with a realistic example near the top of 
the barrier, that it is indeed the choice of a positive sign in the denominator which improves 
the WKB width estimate. The choice of a positive sign also brings (fTOj) in agreement with 

We now move on to discuss one more approach which is popularly used in literature and 
has off late been often used for the evaluation of the half-lives of superheavy nuclei. 



C. Fission model approach 

An approach where the alpha decay is considered as a very asymmetric fission process was 



introduced about two decades ago by Poenaru and co-workers 



10| . This model, also known 



as the super asymmetric fission model (SAFM) has recently been implemented extensively 
for the evaluation of the WKB widths of superheavy nuclei [l^, [2^ . The decay width of a 
metastable state in the SAFM within the WKB framework is given as, 

Tsafm{E) = uP = ^(^1 + e'^y\ (11) 

where, P is the probability of penetration through the potential barrier, E,^ is the "vibra- 



tional energy", K = /;^^ K{r) dr with K{r) = V(2/i/r) {V{r) - Q') and Q' = E + E^. 
The Gamow factor with K over here differs from the W occurring in the equations so far as 
obtained in 0, 1^, due the energy E of the tunneling particle being replaced by + E,y. 
Replacing the assault frequency u from Eq. ([5]) into the above Eq. (fTTj) . we indeed recover 
an expression similar to that of TTPAiE), with the k{r) replaced here by n{r). The SAFM 
width is, 



2 



^safm{,E) — — 



'■2 dr 



-^-2X?^«^^ (12) 

/ri K{r) J 

where we have approximated ( 1 + e^^)^^ ~ e~^^ for sufficiently large K. In fact, if we 
start with the definition oi Ey = (1/2) hu = (1/2) (2 tt z/) (as defined in the SAFM based 



works lol . 



I20I ) and use Eq. ([5]) for = T \ we obtain a theoretical relation for E^,, 



namely, 




dr 



(13) 



'V 



2/i [X 



(2^ltC-) (l/(r) - B - 



One could have of course inferred the above equation f|T3l) simply by the comparison of 
Eqs. (fTTj) and (fT2|) . Provided the potential is known, for a given energy of the tunneling 
particle, Eq. f|T3|) is a transcendental equation for the vibrational energy Ey. 

Comparing the expressions for the widths, Vtpa-, ^dh and Tprdman, away from the ex- 
tremes of the barrier (taken in the limit of large W and negligible da/dE), one can see that 
indeed. 



TsAFM agrees exactly in form with the above formulae for widths. The only difference lies 
in the replacement of E hj E + E^y as mentioned before. It is interesting to note that even 
though the prefactor in front of the exponential in Ttpa arises due to the normalization of 
the WKB wave function, it agrees exactly with the prefactors in Tdh and Tprdman where it 
arises due to the replacement of the assault frequency as in ([S]). 

Coming back to Eq. ( fT3i) . one notices that for a given potential, V{r), in a particular 
decay problem with a given Q- value, E,y can be determined by resolving Eq. ( fT3l) . However, 
starting from the pioneering works of Poenaru and co-workers until some recent ones, E^ is 
fitted to reproduce the half-lives under consideration. There is no mention in these works 
of the fitted value being consistent with (fT3l) . Such a fitting procedure performed without 
a consistency check with Eq. f[T^ would be somewhat ambiguous and outside the spirit of 
a proper WKB calculation. Without the condition (fT3l) . E,y becomes simply a parameter to 
compensate for the mismatch of the theoretical width with experiment. The interpretation of 
El, as a vibrational energy and its addition to the Q value {E = Q) giving, Q' = Q + E^, is 
then not justified. The above point will be clarified with realistic examples of the calculation 
of half-lives of superheavy nuclei in section IV. In the next section, we briefly describe the 
potentials used for the calculations of the present work. 

III. THE ALPHA NUCLEUS POTENTIAL 



TtPa{E) — Tj:,h{E) — T Frmnan{E) 



fl^ r dr 1"^ 



-2 k(r)dr 

e •^'■2 



(14) 



2 Jrx k(r) . 



With the objective of the present work being a critical examination of the various semi- 
classical methods used for the evaluation of alpha decay half-lives, we perform calculations 
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using different available inputs in literature. The potential in ([T]) is written using a double- 



11| . The folded nuclear 



folding model with realistic nucleon-nucleon interactions as given in 
potential is written as, 

Vn{r) = ^ J dvidv2Pairi) Pdir2)v{ri2 = r + - ri, E) (15) 

where pa and pd are the densities of the alpha and the daughter nucleus in a decay and 
v{ri2,E) is the nucleon-nucleon interaction. |ri2| is the distance between a nucleon in the 
alpha and a nucleon in the daughter nucleus, f (ri2, E) is written using the M3Y nucleon- 
nucleon (NN) interaction as in pjj] as, 

v(.,,,E) = 7999H?tlM) _ ,i3, e.p(-2.5|r..|) ^ ^^^^^^^^ ^^^^ 

4 |i'i2| ^-0 |ri2| 

Joo = -276(1 - 0.005 EJA„). 

The alpha particle density is given using a standard Gaussian form, namely. 

Pair) = 0.4229 ea;p(-0. 7024 r^) (17) 

and the daughter nucleus density is taken to be, 

where po is obtained by normalizing pd{r) to the number of nucleons Ad and the constants 
are given as c = 1.07 ^y^fm and a = 0.54fm. The equation (ITS!) involves a six dimensional 
integral. However, the numerical evaluation becomes simpler if one works in momentum 
space as shown in [11]. The constant A is determined by imposing the Bohr-Sommerfeld 
quantization condition ([2]) using the above potential. The number of nodes are re-ex pre ssed 
as n = (G — I) /2, where G is a global quantum number obtained from fits to data 12, 2l| 
and / is the orbital angular momentum quantum number. We shall perform calculations 
with two possible fitted values of G [2])], namely, 22 and 24. The Coulomb potential is 
obtained using a similar double folding procedure with the matter densities of the alpha 
and the daughter replaced by their respective charge density distributions p'^ and p^. Thus, 
double folding the proton proton coulomb potential, 

|ri2| 



V;(r) = / dr, dv2 p^(ri pS(r2) f- . (19) 
J ri2 
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The charge distributions are taken to have a similar form as the matter distributions, except 
for the fact that they are normahzed to the number of protons in the alpha and the daughter. 

One could further improvise the double folding potential by taking into account the 
density dependence of the NN interaction v{ri2). For example, in [22|] a reasonably good 
description of elastic alpha-nucleus scattering data was obtained by assuming a factorized 
form of the density dependence as follows: 

{)(ri2, p„, pd, E) = C i;(ri2, E) /{p^, E) /{p^, E) , (20) 

where, f{px, E) = 1 — {3 p^x^ , with X being either the a or d. The parameters, C and 
(3 were found to be energy independent and C = 1.3 and (3 = 1.01 fm^ for the range of 
analysed data between alpha particle energies of 100 MeV and 172 MeV. Due to the lack of 
much information, for the case of super heavy nuclei, C was chosen to be unity and (3 = 1.6 
fm^ {19!]. We note here that even if the potential is improvised with the density dependence 
of the NN interaction included, the Bohr-Sommerfeld condition should still be satisfied in 
the WKB framework. As we shall see later, the normalization A in (|T5l) is different from 
unity even for the above interaction. 

In order to test the applicability of the formulae iQ and ffTOj) at the top of the barrier, 
we shall examine the case of the / = 2, ^Be resonance which decays 100% into two alphas. 
The nuclear potential for the a-a case is taken to be |ll|. 

y°"(r) = -122.6225 exp(-0.22r^) MeV. (21) 

Since the aim of the present work is to make a comparative study and not fit parameters to 
match the theoretical widths with the experimental ones, the alpha particle preformation 
probability for all the calculations in this work has been taken to be unity. 

IV. RESULTS AND DISCUSSIONS 

The objective of the calculations performed for the superheavy nuclei is to test the sensi- 
tivity of the results to (i) the implementation of the Bohr-Sommerfeld quantization condition 
([2]) which fixes the strength of the potential A in f|T5l) . (ii) the 'fitted' global quantum num- 
ber G appearing in ([21), (iii) the density dependence of the nucleon-nucleon [DD — NN) 
interaction and finally (iv) to verify if the fitted vibrational energies used in literature are 
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consistent with the theoretical Eq. (fT3l) . Some recent works 19| on superheavy elements in 
literature neglect (i), (ii) and (iv) from above. In what follows, we shall see that the exact 
values of the decay widths of the nuclei considered do depend strongly on the ingredients 
(i), (ii) and (iv) and hence any conclusions drawn in works which neglect these aspects of 
the WKB framework would have to be treated with caution. 



1. Sensitivity of superheavy nuclear half-lives to DD — NN and the Bohr-Sommerfeld condition 

In Tables I and II, the half lives ti/2 = ln(2)/r((5), of some currently discovered super- 
heavy nuclei are shown for two different choices of the global quantum number G of the 
Bohr-Sommerfeld quantization condition. Since all the decays discussed in Tables I and 
II take place at energies E = Q which are away from the extremes of the barrier, we use 
T{Q) = Ttpa{Q) which in turn is the same as evaluating Tdh{Q) or T prdmaniQ) as shown 
in (iMl) . Theoretically, the Q value of the decay is defined as the difference of the masses 
of the parent nucleus and the sum of the masses of the alpha and the daughter nucleus 
{Q = Mparent — — Mrf). Wc shall however use the Q deduced from the measured a- 
particle energies, Ea, by applying a standard recoil correction as suggested by Perlman and 



Rasmussen 



and frequently used in literature. With Zp and Ap being the charge and 



mass numbers respectively of the parent nucleus, 

Q = ^ E^ + (65.3 - 80.0 Z^/^ ) IQ-^ MeV . (22) 
Ap — 4 

The results obtained using Eq. ( |T6|) for the nucleon-nucleon (NN) interaction are labelled 
as "free-NN" in Tables I and II. The density dependent NN interaction {DD — NN) calcula- 
tions use Eq. (!20l) instead of f (ri2) in (fT5l) and are also shown in the tables. One can see that 
the introduction of density dependence in the NN interaction reduces the lifetimes ti/2 by 
an order of magnitude as compared to the free NN results. Neglecting the BS condition (i.e. 
A = 1) however, 'increases' the DD — NNti/2 by two orders of magnitude as compared to 
the proper DD — NN calculation using the BS condition with A around 2. Any conclusions 
based on calculations neglecting the BS condition can hence be quite misleading. In an ideal 
case, when the potential K„(r) is known exactly for a particular system, one would expect 
A in (fT5|) which gets fixed by the BS condition to be unity. One can however see that using 
the 'free NN' interaction, the value of A ranges around 0.6 — 0.7, while for the DD — NN 
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TABLE I: Half lives in ms, for G= 22 



Parent Nucleus 

lllSg Q=8.67 MeV 


A" 

0.644 


ti (ms)" 
17638.6 


2.095 


ti(ms)^ 
2794.5 


ti (ms)'^ 
213197.8 


ti (msY 

2^ ' 

288000 


lllHs Q=9.44 MeV 


0.639 


356.3 


2.080 


56.65 


4077.1 


150 


lllDs Q=11.368 MeV 


0.638 


0.021 


2.072 


0.0034 


0.205 


0.17 


llfRg Q=11.36 MeV 


0.639 


0.0439 


2.076 


0.0071 


0.46 


6.4 


27^112 Q=11.3 MeV 


0.637 


0.117 


2.072 


0.019 


1.253 


0.69 


286114 Q=10.35 MeV 


0.634 


102.07 


2.067 


15.61 


1248.3 


400 


293116 Q=10.67 MeV 


0.627 


57.583 


2.049 


8.701 


679.6 


53 


29^118 Q=11.81 MeV 


0.625 


0.416 


2.042 


0.064 


4.594 


1.8 



Tree NN, ''Density dependent NN, '^Density dependent NN witli A = 1, ''Experimental value 



case it is in the range of 2 — 2.3. It is however interesting to note that the value of A hardly 

depends on the mass or atomic number of the parent nucleus for the considered range of 

nuclei. The parameter C in (|20l) was in fact chosen to be unity due to lack of information. 

One could rather choose C = 2 and C = 2.3 leading to a A close to 1 for the DD — NN 

cases in Tables I and II respectively. 

Comparing the numbers in Tables I and II, one observes that the half lives are reduced 

on increasing the value of the global quantum number from G = 22 to G = 24. With the 

primary objective of the work being a comparison of the different WKB approaches, the 

orbital quantum number I was taken to be zero. Note however, that the introduction of the 

Langer modification, namely, /(/ + 1) ^ (Z + 1/2)2 ^j^g radial one-dimensional WKB 

problem introduces an additional turning point near the origin, even for the / = case [14 1 . 

□ 

This detail has also been missed out in some works [19|. 

In passing, we note that in some of the cases like the decay, "^^^114 —>■ a + "^^"^112 for 
example, the lifetime of the daughter nucleus, ti/2(^^^112) = 0.5 ms, is much smaller than 
that of the parent, ti/2(^^^114) = 160 ms (implying that the daughter in the cluster decays 
before the parent can decay). The application of quantum tunneling (which assumes a 
preformed cluster) to such a problem would be somewhat ambiguous with the daughter 
decaying faster; however, one could also argue that the daughter inside the cluster does 
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TABLE II: Half lives in ms, for G= 24 



Parent Nucleus 

lllSg Q=8.67 MeV 


A"* 

0.720 


ti (ms)" 
10511.8 


A" 

2.337 


ti(ms)'' 
1680.8 


1 1 (ms)'^ 

2 ^ ^ 

213197.8 


ti (msY 
288000 


lllHs Q=9.44 MeV 


0.715 


213.1 


2.320 


34.185 


4077.1 


150 


lllDs Q=11.37 MeV 


0.713 


0.012 


2.312 


0.0021 


0.205 


0.17 


llfRg Q=11.36 MeV 


0.714 


0.027 


2.316 


0.0044 


0.46 


6.4 


27^112 Q=11.3 MeV 


0.712 


0.0702 


2.311 


0.0114 


1.253 


0.69 


286114 Q=10.35 MeV 


0.707 


60.57 


2.302 


9.355 


1248.3 


400 


293116 Q=10.67 MeV 


0.699 


34.102 


2.280 


5.213 


679.6 


53 


29^118 Q=11.81 MeV 


0.697 


0.248 


2.272 


0.038 


4.594 


1.8 



Tree NN, ''Density dependent NN, '^Density dependent NN with A = 1, ''Experimental value 



not decay as fast as the free one and the picture is still valid. Indeed, in recent literature 
the tunneling picture is used without considerations of the lifetimes of the parent and the 
daughter nuclei. 



2. Assault frequencies and vibration energies 

We shall now examine the evaluation of the widths using Tsafm of Eq. ( |T2|) . As men- 
tioned before, the only difference in the evaluation of Tsafm as compared to Ttpa is in 
the replacement of Q by Q + E^. The zero point vibration energies, Ey, are usually taken 



from fits 2J] and are given for the superheavy case [19] as, Ey = 0.1045(5 for even-even, 
Ey = 0.0962Q for odd Z - even N, E^ = O.OQOTQ for even Z - odd N and E^ = 0.0767g 
for odd Z - odd N parent nuclei. With an average Q value for the superheavy nuclei around 
Q = 10, one could say that E^ would be of the order of 1 MeV. Such values are however 
not consistent with Eq. (fT3|) . If for example, we provide the above values of E^, from fits as 
an input for the right hand side of Eq. f|T3|) . the outcome (which in principle must be E^, 
itself) turns out to be a much larger energy. For the eight nuclei considered in tables I and 
II, we evaluated the right hand side of Eq. ( |T3i) providing E^, as an input. To be precise, we 
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TABLE III: Comparison of the vibrational energies and corresponding assault frequencies obtained 
from fitted values, Eq. ([23]) and Eq. (l2iD . 



Parent E^ju i^.,test E^^tpa ^fit = {-^) latest = {^^) i^tpa = {-^^^) 
Nucleus (MeV) (MeV) (MeV) (lO^is'i) (W'^^s-^) (lO^^s-^) 



271 c 
106 '-'5 


0.786 


4.130 


6.043 


0.380 


1.998 


2.923 


275 u 
108 -"'5 


0.856 


4.257 


5.987 


0.414 


2.059 


2.896 




1.031 


4.226 


5.903 


0.499 


2.044 


2.856 


fARa 


0.871 


4.220 


5.917 


0.421 


2.041 


2.863 


277112 


1.025 


4.210 


5.913 


0.496 


2.036 


2.861 


286;^]^4 


1.082 


4.184 


5.949 


0.523 


2.023 


2.878 


293116 


0.968 


4.153 


5.910 


0.468 


2.010 


2.859 


294118 


1.234 


4.127 


5.856 


0.597 


1.996 


2.833 



evaluated, 



u.test 



dr 



-1 



(23) 



using the fitted values of E^, defined as E^ju in the equation above. The Ey^test?> are not 
the same as Eujus (as should have been the case due to Eq. (fT3|) ) and are listed in Table 
III. The calculations were done with the density dependent NN interaction in the nuclear 
potential and with the value of A = 1 to perform a comparison with the works which use 
the SAFM. The vibration energy is in fact related to the assault frequency at the barrier as, 
E,y = (1/2) h {2 7T u). In Table III we also list the assault frequencies corresponding to the 
fitted -Ejy's as used in the SAFM models and to the calculated E,^^test values. For comparison, 
we present the assault frequencies appearing in the widths, Ttpa, ^dh and Tprdman which 
are the same in all the three cases (see Eq. ([5])) and label them as i'tpa- The corresponding 



E,. 



jPA IS given as. 



E,. 



TPA 



2/i 



r2 



dr 



-1 



(24) 



'r^ ^J{2f,/h') (F(r) - Q) 

The above calculation of E^^tpa is performed for a density dependent NN interaction and 
including the BS condition. The assault frequencies, vtpa are of the order of 102^s~^ which 
is more like the standard result expected for alpha particle tunneling 25 1. 
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3. Improved WKB formula at the top of the harrier 



Finally, we discuss the result of the application of the improved formulae ([9]) and (fTOj) 
for the decay taking place with an energy close to the top of the barrier. Such examples are 
indeed difficult to find among the alpha decay of nuclei. A suitable one is the decay of the 
®Be (2^) level at 3.03 MeV above its ground state. The experimental width of this state 
is 1.513 MeV with 100% a decay. Using the analytical nuclear potential fl2Tl) . the folded 
Coulomb potential and the usual formula ([3]) for the WKB decay width Ttpa for regions 
away from the barrier gives a theoretical width of 1.232 MeV for this level. However, with 
the barrier height being 3.27 MeV and the Q value of the decay ^Be a + a being 3.122 
MeV, the use of the standard WKB formula Q is not recommendable. Instead, if we use 
Eq. ffTOl) but with a positive sign in the denominator as explained below ffTOj) . the width 
'^^FydmaniQ) = 1-535 McV and is closer to the experimental value of 1.513 MeV. Using the 
formula of P], T*^^{Q) = 1.53 MeV which is again close to the experimental number as well 
as consistent with ^^FrdmaniQ) ■ If use the expression of {4] as it is in (fTOl) with a minus 
sign in the denominator, the width turns out to be 0.63 MeV which indeed worsens the 
result of 1.232 MeV, obtained with the standard WKB formula and also disagrees with the 
rZiQ) - 1.53 MeV^ Since the autho. had a eho.ee of the s,g. in fl and chose the negative 
sign without any particular argument, we guess that the choice should have rather been the 
opposite and the expression should be read with a T + 2hda/dE in the denominator. 



V. CONCLUSIONS 



To summarize, we first performed a survey of the available WKB decay width formulae in 
literature, which were obtained using different models and approaches. After having noted 
the similarities as well as differences in the various approaches, we applied them to the cal- 
culation of the half lives of super heavy nuclei which form a topic of current interest. The 
motivation to apply for the case of superheavy nuclei was also to emphasize the need for 
performing calculations remaining consistently within the spirit of the WKB approximation. 
Following are the main observations of the present work: 

(i) The decay widths of the super heavy nuclei are sensitive to the input of density depen- 
dence in the nucleon-nucleon interaction of the nucleons in the a particle and the daughter 
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nucleus. Since the half lives can reduce by an order of magnitude as compared to the results 
with a free NN interaction, any conclusions drawn in such works regarding the angular mo- 
mentum, values, become model dependent. 

(ii) Conclusions obtained in some recent fission model based calculations of super heavy nu- 
clei neglect the Bohr-Sommerfeld condition which amounts to discarding the semi classical 
nature inherent to the WKB approach. We find once again that the half lives can change 
by orders of magnitude by neglecting this condition. 

(iii) The results as in (ii) above, often seem to be in agreement with data (see Tables 
I and II). However, in comparing theory with experiment one should not get tempted to 
choose an inconsistent approach for the sake of obtaining agreement as is often done in the 
SAFM calculations. 

(iv) The assault frequencies appearing in the fission model calculations are shown to dis- 
agree with three other approaches existing in literature. These frequencies and hence the 
vibrational energies which are fitted in the fission model based calculations are in principle 
inconsistent with the formulae obtained from the standard WKB method. 

(v) Improved formulae for the decay at the top of the barrier are compared and applied to 
a realistic example. We suggest the fiip of a sign in the denominator of the expression f|TO]) 
obtained in the work of Drukarev, Froman and Froman Such a sign fiip brings the results 
in agreement with experiments as well as with Eq. ([9]) for the width from another work . 
The sign fiip is consistent with the theory in ^, since at some point in their derivation one 
encounters a choice of signs. 

To compare theoretical WKB widths with experiment, it is mandatory to perform a 
consistent calculation, taking into account carefully the details like the Langer modification, 
Bohr-Sommerfeld condition and a theoretically derived vibrational energy. We conclude by 
mentioning that any attempt to extract physical information from fitted parameters (such as 
the alpha cluster preformation probabilities or the unknown angular momenta of superheavy 
nuclei) while calculating the half lives within the WKB approximation should bear in mind 
the limitations introduced in the model due to the sensitivities mentioned above. 
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